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ABSTRACT 

All-sky maps of the cosmic microwave background temperature fluctuations are usually 
represented by a spherical harmonic decomposition, involving modes labeled by their 
degree I and order m (where —l^m^ +1). The zonal modes (i.e those with m = 0) 
are of particular interest because they vary only with galactic latitude; any anomalous 
behavior in them might therefore be an indication of erroneous foreground substrac- 
tion. We perform a simple statistical analysis of the modes with low I for sky maps 
derived via different cleaning procedures from the Wilkinson Microwave Anisotropy 
Probe (WMAP), and show that the zonal modes provide a useful diagnostic of possible 
systematics. 

Key words: cosmology: cosmic microwave background — cosmology: observations 
— methods: data analysis 



1 INTRODUCTION 

Observations of the temperature anisotropies of the cos- 
mic microwave background (CMB), particularly those from 
the Wilkinson Mic r owave Anisotropy Probe (WMAP) 
jBennett et al.ll2003l : Iffinshaw et al"1l2009l) . form the foun- 
dations of the re markably su ccessful "concordance" cos- 
mological model (I Coles! 2005). An essential ingredient of 
this model is the assumption that the primordial den- 
sity fluctuations that seeded the formation of galaxies and 
large-scale structure were statistically homogeneous and 
Gaussian. Versions of the inflation scenario based on the 
idea of a single slow-rolling scalar field predict levels of 
non-Gaussianity too small to be observed. On the other 
hand, multi-field inflation models, and models with a non- 
standard kinetic term for the inflaton, may yield larger non- 
Gaussian effects which could in principl e be detected in cur- 
rent or next-genera t ion ob s ervations jBartolo et alj 20021: 
Bemardeau fc Uzanl l2002t iLvth et al] 120031; iDvali etafl 
20041: lArkami-Hamed et alj |2004 lAlishahiha et alj 12004 



Bartolo et al. 



2004; 



...Chen et alj|2007l : iBattefeld fc Battefeldl 
20071 : I Kovama et al. I 120071 ). Analysis of currently avail- 
able WMAP dat a provide strong limits on the level of 
non-Gaussianity dKomatsu et alj 20031: Spergel et alj 120071 : 
ICreminelli et alj 120071: iHikage et al.1 120081 ). On the other 



hand, lYadav fc Wandelt] (|2008l ) recently reported a detec- 
tion of primordial non-Gaussianity at greater than 99.5% 
significance. Further detailed analyses of non-Gaussianity 
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are clearly necessary in order to reconcile and understand 
the various constraints and claimed detections. 



The greatest barrier to the detection of non- 
Gaussianity, or other departures from the framework of 
the concordance cosmological model, is the presence of 
residual foreground contamination or other systematic er- 
rors. Since our own Galaxy emits at microwave frequen- 
cies, the emission from local foregrounds must be care- 
fully cleaned before a map can be obtained that is suit- 
able for analysis. One way of avoiding this problem is to 
cut out regions near the Galactic plane where contam- 
ination is particularly severe, but this throws away the 
advantage of having full coverage. It is therefore impor- 
tant to produce maps that are as clean as possible over 
the whole sky for many purposes. However, such clean- 
ing i s inevitably approx i mate and biases are bound to oc- 
cur jEriksen et aL 20051: iNaselskv. Verkhodanov fc Nielsenl 
|200S| ; IChiang, Naselskv fc Coles! 120091 ). Circumstantial evi- 
dence exists that may be interpreted as being due to the 
presence of residual Galactic foregrounds in the WMAP 
data, or some other artefact of the cleaning process 
l|Chiang, Naselskv fc Colesll2007l ; IChiang et alj|2007l ), but in 
the absence of a more complete characterization of the galac- 
tic emission the situation remains unclear. 



In this Article we propose and test a simple diagnostic 
analysis that offers the possibility of identifying foreground- 
related biases and systematics in all-sky maps of the CMB. 



2 Short & Coles 



2 ZONAL MODES OF CMB MAPS 

The statistical variation of the CMB temperature, T(0,tp), 
over the celestial sphere can be conveniently decomposed 
into spherical harmonic modes: 



1=0 m=-l 



(1) 



where the Yi m (8, ip) are spherical harmonic functions, de- 
fined in terms of the Legendre polynomials, Pj m , using 



?,¥>) = (-!)' 



(2Z + l)(Z-m)! 
4n(l + mj\ 



P im (cos6 l )e 1 



(2) 



and the aim are complex coefficients which can be ex- 
pressed with aim = \aim\ exp(^$ im ) where §i m are the 
phases ilColes et all 12004 IChiang. Naselskv fc Coles! |2004 



12001 120091 ; IStannard fc Colej|2005h . We have adopted the 
Condon-Shortly phase definition in the spherical harmonic 
decomposition. 

The spherical harmonics functions, Yim, can be visual- 
ized by considering their nodal lines, i.e. the set of points 
(9, ip) on the sphere where the spherical harmonics vanish, 
i.e. where Yi m (9, <p) = 0. Nodal lines are circles, which are 
either in the latitude or longitude direction with respect to 
the coordinate system being used; in the case of CMB maps 
this is usually the galactic coordinate system. The number of 
nodal lines of each type is determined by the number of zeros 
of Yi m in the latitudinal and longitudinal directions indepen- 
dently. The associated Legendre functions Pi m possess I— \m\ 
zeros in the latitude direction, whereas the trigonometric sin 
and cos functions possess 2\m\ zeros in the longitude direc- 
tion. Two specific orders m are of particular interest at a 
given I in the context of this work: these denote the zonal 
modes, with m — 0, and the sectoral modes, with m — I. 
In the former case there are no zero-crossings in the longi- 
tude direction, so contours of equal temperature run parallel 
to latitude lines; in the latter the contours run parallel to 
longitude lines. Examples are shown in Figure 1. In the in- 
termediate cases with < m < I there are zero-crossings 
in both directions, giving rise to a patchwork appearance; 
these are usually called tesseral modes. 

Statistically isotropic Gaussian random CMB temper- 
ature fluctuations on a sphere, of the type that result from 
the simplest versions of the inflationary paradigm, possess 
spherical harmonic coefficients (ai m ) whose real and imag- 
inary parts are mutually independent and both Gaussian 
l|Bond fc Efstathioull 19871 : IColes et al.|[2004) . The statistical 
properties of the fluctuations are then completely specified 
by the angular power spectrum, C; , where 



(ai m a*> m <) = Ci 6 U ,6„ 



(3) 



Since T is always real, the complex vectors of the a; m on 
the Argand plane for m < are mirror images of those with 
m > with respect to x axis for even m, and with respect 
to y axis for odd m. Writing 



aim = Xl m + iyir, 



(4) 



we note that the variances of the real and imaginary parts 
of aim for m > are equal 



cr 2 {xim) = o- 2 (yim) = erf = -Cj, 






(5) 
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Figure 1. Illustrative examples of zonal and sectoral modes for 
I = 10; the first example is the zonal mode with m = and the 
second is the sectoral mode with m = 10. 



which depends only on I. The distributions of xi m and 
yim are independent Gaussians with these variances and 
zero mean. The amplitudes \ai m \ for these modes there- 
fore have a Rayleigh distribution with rando m phases 
jBond fc Efstathioulll987l ; IStannard fc Colej2005f ). For m = 
0, the imaginary part of a; m must be zero, so this mode al- 
ways has zero phase $! m . This is because the phase relates 
to the variation around the polar axis only. However, the 
distribution xi m for m = should be equal to that of the 
other real parts at a given I, namely a Gaussian with zero 
mean and stated variance. 



3 STATISTICAL ANALYSIS 
3.1 Preamble 

Statistical analysis of the spherical harmonics of CMB maps 
usually involves using all modes in an equivalent manner. 
However, since the zonal modes in particular have such a 
special relationship to the galactic coordinate system, it 
is worth looking at their properties independently to see 
whether they hold any clues to possible residual contamina- 
tion aligned with the Galactic plane. 

In order to construct a test which involved the smallest 
possible number of assumptions, and in particular avoided 
the need to make estimates of the power-spectrum Ci along 
the way, we focused on the modes with maximum or min- 
imum amplitude at a given I. If all modes are statistically 
equivalent then the different orders m at a given I are equally 
likely to furnish the maximum (or minimum) amplitude 
|aj m |. A preference for modes with m = to display the 
maximum (or minimum) amplitude might therefore be plau- 
sibly interpreted as evidence that the zonal modes are either 
contaminated with residual foreground, or that foregrounds 
have been excessively subtracted. Both of these possibilities 
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are supported by other evidence ([Chiang. Naselskv fc Coles 
20091 ; IChiang et ai1l2007l ; iNaselskv, Verkhodanov fc Nielsen 

200a ) . 



3.2 Data 

The maps analyze d in th is paper are th e 1-year 
jBennett et al.l l2003h . 3 -year (jjarosik et all 120071 ). and 5- 
year ( Hinshaw et al.l I2009T I Internal Linear Combination 
(ILC) maps from the WMA P team, the 1 and 3 year maps 
from Tegmark et al.l |2003l ) and the harmonic ILC map 
by iKim et all (|2008h (hereafter ILC1, ILC3, ILC5, TOH1, 
TOH3 and HILC). We apply three tests to these maps, 
which are described in detail below. 



3.3 Extremal Mode Counts 

For our first test, we analyze the value of m associated with 
the largest (or smallest) amplitude |a; m | for a given I in the 
ranges [0, 10] and [0, 20]. Since the a; m for all orders m at a 
given I should have the same variance in the null hypothesis 
(of a stationary Gaussian random field), the minimum or 
maximum value of \ai m \ should not occur preferentially at 
any particular value of m. We use a very straightforward 
method to establish whether this is the case. For each map 
we simply count the number of occurrences (i.e. numbers 
of separate degrees I within the range analyzed) for which 
the minimum or maximum value amplitude occurs at m = 
(the zonal mode) or m — I (the sectoral mode) . These counts 
are recorded in Table [1] What we have done is to identify, 
at each degree I, the value of m which has the minimum (or 
maximum) value of |o; m |. If the minimum is at m = then 
this contributes to the count in Column 3 of the table; if the 
maximum is m = it contributes to the count in Column 4; 
likewise occurrences of extrema at m — I contribute to the 
counts in the following two columns. For results where the 
maximum I — 10 the counts are out of a total of 11 and for 
maximum I = 20 they are out of a total of 21. 

We have restricted our analysis to low I modes partly 
to keep the computational cost of the simulations down but 
mainly because the high I modes are known not to be clean 
anyway so it would not reveal anything interesting to find 
zonal or sectoral anomalies among them. We have also given 
results for the sectoral modes for comparison, but there are 
no significant anomalies associated with them and we shall 
not discuss them further in this paper. 



3.4 Significance Levels 

It is not a trivial matter to calculate significance levels ana- 
lytically for this test because the number of available orders 
m increases with I. It is obviously much more probable that 
the minimum amplitude is at m = for I — 2 than for I — 20 
under the null hypothesis. Assessing the significance of the 
number of occurrences of zonal or sectoral extrema involves 
a messy exercise in combinatorics. However, we can finesse 
this difficulty by instead comparing the actual maps with 
simulations constructed according to the Gaussian assump- 
tion described in Section 2. 

It is also possible that the cleaning process used to re- 
move foreground contamination from the raw observations 
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Table 1. This table shows, for the various cleaned all-sky maps 
described in the text, the maximum value of the degree I consid- 
ered, and, in the following four columns, the number of times (i.e. 
number of values of I) for which the minimum (or maximum) am- 
plitude ai m is at m = (zonal mode) or m = I (sectoral mode). 
For example, the m m i n = and max I = 20 result for TOH3 
means that 11 out of 21 different values of the degree I £ [0,20] 
have the lowest value of the amplitude |a; m | at m = 0. 



might in any case introduce some sort of bias into the sta- 
tistical distribution of amplitudes or induce correlations be- 
tween different modes. To circumvent this difficulty, as well 
as the one noted in the previous paragraph, we therefore 
base the confi dence levels on our te st on a set of simulations 
performed bv lEriksen et al.l (|2005l ). which takes "raw" sky 
maps, generated assuming Gaussian fluctuations, adds sim- 
ulated foregrounds and then recovers the signal using the 
ILC methodology. The simulated maps we use therefore al- 
ready take into account any "artificial" correlations that the 
ILC process may generate. These are particularly useful as 
they allow us to assess whether any anomalies we do actu- 
ally find must be above the level known to be introduced by 
the ILC cleaning process. 

In order to calculate significance levels of the results in 
Table Q] we used an ensemble of N = 10000 independent 
Monte Carlo realizations of Gaussian skies, contaminated 
with foreground and then cleaned according to the ILC pre- 
scription as described above. We use these simulations to 
construct empirical distributions of the count statistics dis- 
played in the previous table and from these we compute the 
empirical significance levels shown in Table [2] 

The corresponding probabilities of the results in Table[2] 
were calculated using the Monte Carlo simulations detailed 
above by counting the number of occurrences of each con- 
figuration shown in Table [1] Note that since we are using 
N — 10000 independent simulations we expect the results 
to be affected by Poisson fluctuations at the level of order 
y/N, which means we expect the probabilities in Table [2] 
only to be accurate to about 1% or so. 

For the case of zonal maxima, it is clear that there no 
significant results above 2a (i.e. the 95% level), but the zonal 
minima do show a significant result in the TOH3 and HILC 
maps for I < 10. It is also interesting to note that the TOH3 
map gives a higher significance level than the TOH1 map, 
as does the ILC5 map compared to the corresponding 1 yr 
map ILC1. 

It is worth stressing that each column represents a single 
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Table 2. Monte Carlo estimates of the probabilities of the ex- 
tremal mode counts, i.e occurrences of m ox t = or I, for the 
various CMB maps shown in Table [T] These are computed by 
forming the empirical distribution of the counts over an ensem- 
ble of simulated skies and counting what fraction of the ensemble 
gives the results obtained for the real maps. For example, in the 
case of the m m i n = and max I = 20 result for TOH3, we find 
that, out of 10000 simulations, 9844 have less than 11 (from Table 
[TJ occurrences of minimum amplitudes at m = 0. Given the prob- 
able sampling accuracy of around one percent, we have rounded 
the results. Note that these are discrete distributions, so the zero 
percentages do not necessarily indicate cases of exceptional signif- 
icance, just that it is not possible to have less than the observed 
number of occurrences. In other words, we should treat this as a 
one-sided statistical test. 



statistical test using all the information contained in the 
modes with degree / in the range analyzed, and therefore 
represent genuine a priori significance levels. For example, 
in the case of the m m i n = and max I = 20 result for 
TOH3 described in the caption for Table [5] the probability 
that a genuinely Gaussian CMB sky processed in the way 
we described above, would produce as many zonal minima 
as observed is only (100 - 98.4) = 1.6%. 

3.5 Mode Variances 

We performed a second investigation in order to locate the 
source of the apparent deficit in the zonal modes. The re- 
sults above suggest that the amplitudes when m = are 
low compared to simulations, so this test considers whether 
this is because the variance of the amplitudes is lower for 
m — than for all other m. Since there is no imaginary part 
to the amplitude when m = 0, just the real parts of the aim 
(i.e. the Xi m ) are considered. The variances of the xio for the 
ILC5 and TOH3 maps are calculated and the correspond- 
ing probabilities are calculated by comparing these values 
to the variances of the Xio from the simulations. The results 
in Table [3] provide evidence that the variance of the xio in 
the ILC5 and TOH3 maps is significantly smaller than in the 
simulations. For comparison, the same calculation is also ap- 
plied 'for all m' and 'for all m not equal to 0'. The results 
for all m show that the variance in the TOH3 and ILC5 
maps is lower than in the simulations, which is consistent 
with the previously reported low variance in the combined 
Q + V + W map bv iMonteserin et alJ (|2008h . This leads to 
the further question of whether the zonal mode amplitudes 
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Table 3. Percentage of simulations, for given m, where the vari- 
ance of the iE; m is greater than the specified map. For example 
take the TOH3 result for m=0. We found that, out of 10000 sim- 
ulations, 9757 had a variance of the x;o greater that the variance 
of the xio in the TOH3 map. Since the Poisson fluctuations are 
of order 1% we have rounded the resulting percentage. 




Figure 2. Reconstructed CMB map for WMAP ILC 5 yr map 
using only the m = modes for I = [0, 20] i.e. this map is the 
WMAP ILC 5 yr map which has had all but the |a;o| set to zero. 
The colour scale is marked in milliKelvin. 



are entirely responsible for this low variance. Table [3] shows 
that by removing the m — amplitudes the number of sim- 
ulated maps with a variance greater than the test map falls 
dramatically, suggesting that indeed the low variance of the 
m = amplitudes does has a notable affect on the overall 
variance of the amplitudes. 



3.6 Map Structure 

It is interesting to reconstruct what the CMB sky would 
look like if it only contained the zonal modes, as these are 
the ones that appear to have anomalous properties; the re- 
sult in Figure [2] shows that the deepest minima are just 
below and slightly above the Galactic plane. For our third 
test we considered whether these minima (and also the max- 
ima of the map) are of abnormal amplitude compared to the 
simulated maps; this test is done in pixel-space rather than 
using the a; m . Because the resolution of the maps is quite 
low, we restricted our analysis to the parts of the map that 
are well-sampled (i.e. we neglected the Galactic pole areas). 
Monte Carlo simulations were run to estimate the proba- 
bility of a maximum or minimum as seen in the map. The 
resulting probabili t ies ar e shown in Table [4] Similarly to 
lLarson fc Wandeltl l|2004h , who find that the hot and cold 
spots of the separate WMAP frequency maps are not hot 
and cold enough, we find that the maxima (and minima) 
of the zonal maps are not as high (and low) as expected 
compared to the simulations. The trend reinforced again in 
this table seems to be that the anomalous result becomes 
increasingly significant as the maps supposedly improve. 
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Table 4. Probability of the observed maxima and minima in 
zonal maps, as derived from Monte Carlo simulations. For ex- 
ample take the results for ILC5: out of 10000 simulations, 8124 
have a minimum 'temperature' in the m = only map which is 
less than that observed in the ILC5 m = map (see Figure [2L 
Similarly, out of 10000 simulations, 9902 have a maximum 'tem- 
perature' in the m = map which is greater than that observed 
in the ILC5 m = map. 



4 DISCUSSION 

We have presented a simple statistical analysis based on 
properties of the zonal modes of cosmic microwave back- 
ground maps, i.e. those aligned parallel to the Galactic 
plane. An application of the test to various cleaned CMB 
maps gives interesting results. At the 95 per cent level, 
no significant anoma li es appear in the WMAP ILC maps 
(iBennett et all 120031 ; Ijarosik et ail 120071 ; iHinshaw et~aH 
2009) but there seems to be a significant tendency in some 
other maps |Tegmark et al.ll2003l ; iKim et al.l 120081) to have 
zonal modes with systematically lower amplitudes than 
would be expected in the concordance model. Intriguingly, 
the maps that provide the most significant departures from 
the behavior expected under the null hypothesis are those 
based on later issues of the WMAP data. We have also 
showed that the low variance of the maps is more significant 
where you consider only the zonal amplitudes. This shows 
that zonal modes have a notable contribution towards this 
low variance. Finally, we considered the distribution of max- 
ima and minima in zonal maps which again reinforced the 
earlier finding that the anomalies increase with later data 
releases. The question of whether this is due to increasing 
over subtraction, or whether decreasing noise is revealing an 
underlying anomaly such as a systematic problem, cannot 
be distinguished with current data. 

Of course the maps themselves are not statistically in- 
dependent. Indeed, if the cleaning processes involved were 
perfect then they would all be identical. The different re- 
sults we have found for the different maps are attributable 
to the nature of the zonal modes and their extra sensitivity 
to the structures associated with the Galactic Plane. The ap- 
pearance of significant anomalies in some maps rather than 
others is not a statistical fluke but is clear evidence that 
some cleaning methods leave artifacts in the distribution of 
mode amplitudes. 

It must be noted that the probabilities we quote of 
around 98 to 99 percent are not overwhelming, so the re- 
sults we have obtained are indicative rather than decisive. 
This is not surprising, given the relatively small number 
of modes we used. However, we repeated the analysis for 
a coordinate system aligned with the Ecliptic, rather than 
Galactic plane, and found no significant results at all. This 
lends further credence to our interpretation of the out- 



come of our analysis in terms o f an effect related to over- 
subtraction of Galactic emission dChiang. Naselsky fc Coles! 
120091 : iNaselskv. Verkhodanov fc Nielsen! 12008ft . which be- 
comes increasingly pronounced with each data release. A 
more definitive result will have to wait until more detailed 
foreground subtraction can be attempted, such as will be 
the case with the Planck satellite. 
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